The effect of the dust size distribution on asteroid polarization 

Joseph Masiero 1 , Christine Hartzell 2 , Daniel J. Scheeres 2 

ABSTRACT 

We have developed a theoretical description of how of an asteroid's polarization- 
phase curve will be affected by the removal of the dust from the surface due to a size- 
dependent phenomenon such as radiation pressure-driven escape of levitated particles. 
We test our calculations against new observations of four small (D ~ 1 km) near- 
Earth asteroids [(85236), (142348), (162900) and 2006 SZ 217 ] obtained with the Dual 
Beam Imaging Polarimeter on the University of Hawaii's 2.2 m telescope, as well as 
previous observations of (25143) Itokawa and (433) Eros. We find that the polarization 
of the light reflected from an asteroid is controlled by the mineralogical and chemical 
composition of the surface and is independent of dust particle. The relation between 
the slope of the polarization-phase curve beyond the inversion angle and the albedo of 
an asteroid is thus independent of the surface regolith size distribution and is valid for 
both Main Belt and Near-Earth asteroids. 

Subject headings: Minor planets, asteroids; Polarization 



1. Introduction 



The scattering of light can induce a polarization to an otherwise unpolarized photon flux. In 
the case of asteroids, the dust and regolith on the surface act to instill a linear polarization on the 
reflected sunlight. The d egree of this polarizatio n changes with phase angle, albedo, spectral type, 
and surface composition (jMuinonen et al.ll2002al ) . While single-particle scattering will always result 
in polarization perpendicular to the scattering plane ("positive" polarization), polarization in the 
scattering plan e ("negative" polarizatio n) is o bserved for very s mall phase angles for a wide range 
of objects (e.g. iDollfus Zellnerlll979l . etc.). IShkuratov et al.l (|1994h present a review of possible 
physical causes of this effect, the most promising of which are coherent backscattering models (e.g. 
Muinonenlll989l ) which also can describe the photometric opposition effect. 



In this paper, we determine the effects of different size distributions of surface dust on the 
observed relationship between the polarization of the reflected light and the phase angle of the 
asteroid (that is, the Sun-asteroid-Earth angle). A variety of processes can cause dust loss on 
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asteroids. Electrostatic charging of surface dust particles by Solar radiation, for example, can 



levitate smaller particles meters above the surface (jColwell et al.ll2005l . and reference therein). For 
small bodies these effects will result in the preferential loss of the smallest dust grains, effectively 
truncating the original dust size distribution. Here we test six cases of small near-Earth asteroids 
(NEAs), for which the dust loss should be more pronounced, and measure the effect of dust loss on 
the polarization-phase relation. 



2. Grain size limits 

The surfaces of most atmosphereless bodies in the inner Solar system are composed of a loose 
regolith ranging in size from dust t o boulders. On the moon the smallest particles are dust down to 



0.02 jum in size (jTavlor et al.ll2009l ) while for smaller asteroids such as (25143) Itoka wa the regolith 



appe ars significantly coarser, with an apparent lack of anything smaller than 1 mm (jFuiiwara et al 



2006). With a surface gravity significantly higher than all but the few largest asteroids, the dust 
distribution on the moon sets a strong lower-limit on the dust particle sizes we will consider in this 
study. As we will show in § [3J, the polarization-phase relation is insensitive to the upper-limit of 
the grain size distribution as long as that limit is > 10 (im. 

Dust particles lofted from the surface of an asteroid by micro-meteoroid bombardment or 
electrostatic levitation are immediately subject to Solar radiation pressure effects. This pressure 
can put the particle into orbit around the parent body or allow the particle to escape into a 
heliocentric orbit. A lofted particle's motion can be described by Eq [I] which includes the effects of 
both gravity due to the asteroid (where G is the universal gravitational constant, M the mass of the 
asteroid, and r the position vector of the particle in x, y, z space) and Solar radiation pressure (g s ). 
The radiation pressure is assumed to act in the anti-Sun direction, x, with no time dependency 
(i.e. the x-axis corresponds to the Sun-asteroid line and the orbit of the asteroid around the Sun 
is not included in this approximation). The Solar radiation pressure is described by Eq[2j where 
p v is the albedo of the particle (assumed to be zero initially), C s is a constant of Solar radiation 
(1 x 10 7 g km 3 /(s 2 cm 2 )), B is the mass-to-area ratio of the particle, and R is the distance to the 
Sun. The motion of the particle can then be described by the potential given in EqO 



M 

r = -G—r + g s ii (1) 

9s = BR2 (2) 
M 

U = G 1- <? s (x ■ r) (3) 



A mathematical analysis of th e orbital dynamics of test particles at an asteroid has been carried 
out previously by several authors. iDankowiczl (| 19941 ) developed a precise limit on semi-major axis 
of an orbiting particle under the assumption that the asteroid was stationary. IScheeres &; Marzari 
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(|2002l ) incorporated the motion of the asteroid about the Sun and developed a sufficient condition on 
semi-major axis for a par ticle (in this case a spacecraft) to be trapped in orbit about a small body. 
Bvram &: Scheeresl (J2009) combined elements of these two earlier studies to derive more precise 
limits on semi-major axis of an orbiting test parti cle accounting for its motion. For purposes of 
this paper the analysis in iBvram fc Scheeresl (|2009l ) is too detailed to be summarized simply, and 
thus we will only restate the results of these earlier papers. 



From the analysis by iDankowica (| 19941 ) we find that a simple limit exists for the semi-major 
axis of a lofted particle that generally ensures escape from the small body if exceeded. Mapped 
into our notation we find that escape will occur if the semi-major axis is larger than a escape : 



''escape 




(4) 
(5) 



assuming an albedo of p v = 0. Using a very different approach, IScheeres Marzaril (|2002l ) devel- 
oped a sufficient condition for stability about a small body and found that escape is not immediately 
possible if the semi-major axis is less than o-min'' 



(6) 



In the general problem it is possible for the semi-major axis to evolve, but we ignore this effect 
for this derivation. We note that these limiting semi-major axes scale linearly with the small body 
distance from the Sun. Thus, escape is most likely to occur at perihelion. 

In order to calculate the mass-to-area ratio (B) of a particle, we assume that it is spherical 
with a circular reflective surface area, giving B = | pd dd, where pd is the dust density and dd is the 
diameter of the dust particle. Note that these assumptions are conservative as spheres maximize 
the volume-to-area ratio, and thus these calculations provide a lower limit to the cross-sectional 
surface area of particles escaping. Additionally, we are ignoring any timescale for particle loss once 
in orbit, assuming that sufficient time has passed to clear all particles in the perihelion size- limit. 
From this we can expand Eq[2]to find: 



3{l+p v )C s 
2 p d d d R 2 



(7) 



Substituting Eq[7]into Eq[4]and assuming that the particle of interest rests on the asteroid's 
surface (2a = D ast /2, where D ast is the asteroid's diameter), we can develop an expression that 
relates the minimum particle size expected to be found on an asteroid's surface to the physical 
characteristics of the asteroid: 



< 



3C S 



Pv + 1 



TrG D ast p as tp d R 2 



(8) 
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Equation [8] shows that as the size of the central body decreases or as perihelion distance 
becomes smaller, the minimum surface particle size increases. This can be seen in Fig [H which 
shows with thick lines the relationship given in Eq [8] (normalized by the central body diameter) 
assuming the central body and the dust particle have the same density and an albedo of p v = 0.2. 
We note that dust should have densities up to a factor of 2 more than the bulk asteroid density, 
however this is mitigated by the assumption of spherical shapes for particle size. For asteroids 
spinning near their disruption rate we can assume that o = D as t/2, resulting in a minimum particle 
size that is a factor of four larger than predicted by EqEl 

Table [1] shows the particle size limit for (433) Eros and Itokawa, ignoring any rotation. Note 
that the minimum particle size is inversely related to the asteroid density, diameter and distance 
from the Sun. If we consider the migration of Itokawa from the Main Belt to its current location, 
we see that as the asteroid approaches the Sun, increasingly larger particles are able to escape. If 
the polarization of light reflected from an asteroid is dependent on the size of the particles on the 
surface we would expect NEAs to show changes in their polarization-phase curve as a function of 
the predicted particle size. 



Polarization Effects 



The coherent backscatter mechanism, a second order light-scattering effect, has been shown to 
be the dominant cause of both the observed polarization signatu re of light reflected from asteroids 
at small phase angles as well as the photometric opposition effect (jMuinonen et al.ll2002al ). As phase 
angle increases this transitions to a regime dominated by first-order single-particle s cattering. Using 



numeri cal techniques to calculate the scattering of individual surface elements, iMuinonen et al 
(|2002bl ) find that the polarization of the light can be approximated as: 



P a 
In 



n 



1 



(k d a) 



0) 



n + lj 2[l + (kda) 2 } 

where P is the percent polarization, a is the phase angle in radians, n is the index of refraction, 
k = 2? is the wave number, and d is the scatter separ ation distance. This assumes a single spacing 
for all scattering particles. Other recent models (e.g. iBoehnhardt et al.l 12004 ; iBagnulo et al.l l2006) 
have assumed two discrete scattering components co-mixed at different weightings to fit observed 
polarization phase curves. 

In this work we will instead consider a continuous distribution of grain spacings. The distri- 
bution of voids in a si licate glass follows an approximate power-law falloff for sizes beyond some 
critical peak void size (jMalavasi. et al.ll2006l ). We begin by assuming that the spacing between dust 
on an asteroid's surface also follows a power law, and that the voids have the same minimum and 
maximum characteristic sizes as the dust. This implies that chemical and mineralogical effects are 
fully incorporated into the size of the constituent grains, which we discuss further in § [6j 



Bottke et al. 



(2005) show that for their model NEA population the size distribution is fit by 
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the power law N oc D for sizes from 100 m to < 10 cm. Using a large-area camera network, 



Hallidav et al.l (|1996l ) find that the mass distribution of fireballs (0.1 < M < 12 kg) entering 
the Earth's atmosphere has an exponential power between —0.5 and —1, translating to a size 
distribution power between —1.5 to —3. We focus on this range of size distribution powers for our 
simulations. 

By integrating Eq[9]over all grain sizes from d m i n to d max we obtain: 

p _ A ^£_(wy_j*j £ £_, iH (10) 



2n \n + lj 2 [1 + (k da) 2 } 



where A is a variable used to fit the depth of the negative polarization, p is the exponential power 
of the grain-size distribution, and d m i n and d max are minimum and maximum sizes of the grain-size 
distribution in microns, respectively. 

If we consider only a single wavelength (in this case A ~ 600 nm, meaning k = 10 ^m -1 ) we 
find that the arbitrary variable A reduces to a simple expression of c^nm? giving: 

P = 1.2 d^n / — " ( Z—^ ) n m i —c j „^ dP dd (11) 

" drain 



2n \n + lj 2 [1 + (10 d a) 2 ] 



Figure [2^, shows the change in the polarization-phase relation for fixed n, d m i n and p. As 
d max increases its effect on the curve diminishes, meaning that the polarization is insensitive to the 
largest particle size for general distributions. For all other simulations we fix the maximum particle 
size at 10 jiva. to limit computation time. In Fig [2)3 we plot the effects of a variable d m i n . As a 
normalizing agent, d m i n has a strong effect on the slope of the relation beyond the inversion angle 
(the point at which the polarization-phase curve recovers to zero), but also on the shape and depth 
of the negative branch of polarization. 

A variable index of refraction most significantly affects the location of the inversion angle ao- 
As shown in Fig[2h small changes in the index of refraction result in major changes to the inversion 
angle and depth of negative polarization, with only minor effects on the slope or shape of the 
negative branch. Finally, the effect of variations in the power of the size distribution are shown in 
Fig[2]l. For the ranges of powers considered there are only minor changes to the slope and inversion 
angle, and only in the extreme case is there a significant change to the negative polarization branch. 

From Eq [11] we can calculate the properties of the polarization-phase relation in terms of 
physical constants only. We find that the slope for a > ao is described by: 

_ -1-275 , p+1 _ p+1 , 

n(p+ I) ' 

where h is the slope (in percent polarization per degree). As d max becomes large (>> 10 fim) it's 
effect on h becomes negligible. 
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Similarly, for n > 1.5 and d max > 10 /im we can describe the location of the inversion angle 
approximately as: 




1 



a « a n 



n + lj 2 (10 dr, 



(13) 



where ao is the inversion angle in radians. 



Cellino et al.l (119991 ) derive updated constants for the relationship between polarization slope 
h and geometric albedo p v : 



l ogPv = -1.118 log h- 1.779 



(14) 



Using Eq[l2]and assuming that d max >> 10 /jm we can now describe the geometric albedo of an 
asteroid as: 



logp v = 1.118 log (n {-p - 1) d min -P~ 2 ) - 1.896 



(15) 



4. Observations 



In order to test the effect of dust depletion on the polarization of small NEAs, we conducted 
observations of four targets over a five week period at the end of 2008 using the Dual-Bea m Imaging 
Polarimeter (DBIP) on the University of Hawaii's 2.2 m telescope (jMasiero et al.l 120071 ). DBIP is 
designed to measure both linear and circular polariz ations of asteroids a nd other point sources at 
magnitudes 10 < V < 17 with errors less than 0.1% (jMasiero et al,ll2008l ). 



Targets were chosen from all small NEAs (D ~ 1 km) that would be brighter than V = 18 
over the observing window. Four asteroids had polarimetric measurements of sufficient quality to 
compare with the theory discussed above. Asteroids (85236), (142348), (162900) and 2006 SZ 2 i 7 
were measured at a range of phase angles sufficient to determine both qo and h. Table [2] gives the 
target name, UT date of observation, apparent V magnitude, exposure time, number of exposures, 
Solar phase angle a, measured linear polarization, and angle of polarization referenced to the vector 
orthogonal to the Sun-object-Earth scattering plane. 

In Table [3] we present for each asteroid the fitted inversion angle and slope (considering only 
phase angles a > 15°) as well as perihelion distance (q), absolute magnitude (H„), albedo (as 
calculated from Eq [T4"|) . diameter (D), n and d m i n calculated from Eqs [T2l and [T3l (assuming p = —3 
and d max = 10), and the range of predicted minimum particle sizes remaining on the surface (d prec [) 
depending on the assumed density (between 1.5 and 3.0 g cm -3 ) follow ing § [H Included i n the 



table are the same values for Eros (jZellner fc Gradidll976l ) and Itokawa (|Cellino et al.ll2005l ). two 
NEAs that have been extensively studied bot h remotely as well as in situ by spacecraft visits. 
A measured density for each of these objects ([Cheng] 120021 : lAbe et al.ll2004l ) allows us to make a 
singular prediction for minimum surface particle size. The d pre d range or value for each of the six 
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asteroids is shown in Fig [TJ Scaling the levitation calculations from IColwell et al.l (|2005l ) to smaller 
asteroids we find our predicted minimum particle sizes are within the range of particle sizes that 
are capable of being levitated, assuming the particle's surface electric potential is equal to that of 
the asteroid's surface. 

Figure [3] shows our observations and fits for each asteroid, as well as for the literature data 
for Itokawa and Eros. Note that the fits for 2006 SZ217 and (85236) are under-constrained, thus 
values for the slope have large errors and the inversion angles are only approximate. In all cases, 
the polarization-phase relation beyond «o was assumed to be strictly linear. 



5. Results and Discussion 

Calculations of the dust retention for the NEAs investigated predict over two orders of mag- 
nitude in variation of the minimum particle size. We find no correlation in our data between the 
dmin calculated from our measured slope and the predicted minimum particle size. If the indi- 
vidual dust particles were unique scattering elements we would expect the opposite to be true. 
Thus the physical size of the dust cannot be the primary cause of polarization or even a significant 
contributor. 

In fact, (85236), (162900), 2006 SZ217, Itokawa and Eros all show slopes and inversion angles 
consistent with those measured for large S-type Main Belt asteroids (MB As), while (142348 ') shows 



polarization similar to what is seen for large C-type MBAs (see, e.g lMuinonen et al.ll2002al ). From 
this we infer that the polarization of light reflected from an asteroid is constrained primarily by 
the way the local mineralogy and chemistry dictate the formation of spaces in the mineral matrix. 
The calculated value of d m i n then represents not the size of the smallest particle on the surface, but 
instead a characteristic minimum size of mineralogical voids. For asteroids with similar spectral 
signatures, and thus similar surface chemistries, we expect to see correspondingly similar polariza- 
tion properties. Asteroid polarization-phase curves should show no variation with changing regolith 
conditions for asteroids of the same composition. 

The polarization-albedo relation (Eq [T4"|) then does not depend on surface regolith size dis- 
tribution and so is identical for NEAs and MBAs. Variations in surface regolith and ob serving 



geometries between NEAs a nd MBAs require thermal models tailored for each population ([Harris 



1998 



Wolters Greenll2009l ) which in turn affects infrared albedo determinations. Polarimetrically 
determined albedos are therefore immune to this complication and can be robustly applied to any 
asteroid regardless of location in the Solar system. 
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6. Conclusions 

Through modeling of the effects of a distribution of spacings between scattering elements on 
the observed polarization-phase relation, we have shown that the relation is independent of the size 
of the particles on the surface. Instead, the slope of the relation beyond the inversion angle (and 
thus the spacing parameter d) remains constant for a range of objects with dust-retention sizes 
spanning two orders of magnitude. The polarization-phase relation is therefore consistent between 
MBAs and NEAs of the same spectral type. An asteroid's surface chemistry and mineralogy are 
the dominant cause of the polarization-phase relation, and this relation is equally applicable for 
both NEAs and MBAs. 
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Table 1. Calculated minimum particle sizes from the dynamical model 





Eros 


Itokawa 


Itokawa 


Itokawa 




Perihelion 


Main Belt 


Aphelion 


Perihelion 


Asteroid Diameter (km) 


16.84 


0.33 


0.33 


0.33 


Distance from Sun (AU) 


1.133 


2.500 a 


1.695 


0.953 


Density (g/cm 3 ) 


2.7 


1.9 


1.9 


1.9 


Minimum Particle Size (/im) 


0.5 


10.3 


22.4 


71.0 



Diameters and orbital distances from JPL/Horizons; density values from Ichcng (2002); Abe et al. (2001); a assumed value 



Table 2. Asteroid Observations 



Asteroid 


UT Obs Date 


V mag 


T e xp (s6c) 




Q 


Linear % Pol a 


P 


2006 SZ217 


2008-11-18 


15.6 


300 


18 


19.8° 


0.1 ± 0.2 


9± 10 




2008-12-3 


16.0 


240 


12 


21.7° 


0.3 ±0.1 


11 ± 7 




2008-12-23 


17.4 


300 


12 


36.9° 


1.4 ±0.5 


176 ± 15 


162900 


2008-11-18 


14.9 


300 


12 


5.3° 


-0.5 ±0.1 


86 ± 13 




2008-12-3 


14.9 


120 


12 


16.4° 


-0.2 ±0.1 


90 ±9 




2008-12-23 


15.1 


100 


18 


30.7° 


1.1 ±0.1 


177 ±7 


142348 


2008-11-18 


16.2 


300 


6 


39.3° 


5.1 ±0.5 


177 ±4 




2008-12-3 


16.1 


270 


18 


32.6° 


2.4 ±0.1 


2±2 




2008-12-23 


16.1 


270 


6 


17.1° 


0.1 ±0.1 


20 ± 11 


85236 


2008-12-3 


16.3 


300 


18 


60.3° 


3.8 ±0.1 


179 ±1 




2008-12-23 


17.7 


300 


12 


46.6° 


2.5 ±0.5 


177 ±9 



a quoted errors are la statistical errors; systematic errors are ss 0.05%. 



Table 3. Derived Physical and Dust Properties 



Asteroid 


«o (deg) 


h 


q(AU) 




Pv 


D (km) 




n 


drain (/iHl) 


dp re d (/Xm) 


2006 SZ217 


~ 18 


0.08 ±0.03 


1.200 


17.4 


3 +0 - 2 


0.8 ±0.2 




1.64 


°-2 


8—32 


162900 


19 ±3 


0.09 ±0.01 


1.215 


15.8 


0.24 ±0.03 


1.9 ±0.1 


1.7 


±0.1 


4 9+0-9 

*- Z -0.7 


3—13 


142348 


17 ±2 


0.15 ±0.01 


1.119 


18.2 


0.14 ±0.01 


n 01 +0.03 


1.60 


±0.06 


2.6 ±0.2 


8—32 


85236 


~ 18 


0.09 ±0.04 


0.850 


18.5 


o.2i° ;? 


0.6 ±0.2 




1.65 


4+j 


20—79 


Itokawa 


20 ± 0.2° 


0.084 ± 0.001° 


0.953 


19.2 


0.265 ±0.002 


0.373 ±0.001 c 


1.72 


±0.01 


4 41+0 03 

4 - 41 -0.05 


62.8 


Eros 


21.7 ±0.9 b 


0.107 ±0.002 t 


1.133 


11.2 


0.203 ± 0.004 
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Fig. 1. — Ratio of the minimum dust particle size on the surface to the diameter of the asteroid, as 
related to the asteroid's diameter times the distance to the Sun. The solid and dotted lines show 
the relationship for constant densities of 1.5 and 3.0 g cm -3 , respectively, when assuming an albedo 
of 0.2. The points indicate the positions of NEAs considered in this study with known densities, 
while the ranges represent NEAs modeled with densities from 1.5 to 3.0 g cm -3 . Deviations of the 
calculated points from the lines of constant density are due to slight differences in the measured 
albedos. 
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Fig. 2. — Simulated polarization-phase curves for variations in the value of (a) the maximum 
grain size, (b) the minimum grain size, (c) the index of refraction, and (d) the power of the size 
distribution. The curves show the relation between percent of incident light polarized, referenced 
to the normal to the scattering plane, and the Solar phase angle (a) measured in degrees. The 
dotted curve in each box shows the polarization described by Eq[9]for n = 1.7, (k d) = 20, and a 
normalizing constant of 50: a good approximation for the polarization from a typical S-type Main 
Belt asteroid. 
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Fig. 3. — Observations from Table [21 as well as the best fitting linear polarization-phase relation 
for phase angles a > 15° (dashed line). Fitted slopes and inversion angles are give n in Table [31 
Note that 2006 SZ217 and (85236) are under-constrained. Data for Itokawa is from ICellino et al.l 
(2005) while data for Eros is from lZellner Gradid (119761 ). 



